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FOREWORD 


This document reports on the turbulence measurements performed In a swirl- 
ing, confined jet flowfield. The experiments were supported by Oklahana 
State University (Purchase Requisition 31256) as part of NASA Grant No. 
NAG 3-74, ''Investigations of Flowfields Found in Typical Combustor Geo- 
metries," The experiments conducted in the Dynamics Technology facility 
were performed with a triple hot-wire probe utilizing a specially con- 
structed probe drive fabricated by students at Oklahoma State Univer- 
sity. Included as Appendix A of this report is AIAA paper No. 82-1262 
which reports on turbulence measurements performed in a confined jet at 
Oklahoma State using a six-orientation single hot-wire technique. A por- 
tion of the analysis of the data in that paper was performed at Dynamics 
Technology under the noted subcontract. 

The authors would like to acknowledge the assistance of Professor David 
Lilley of OSU and especially his students in support of our experiments. 
The financial support of NASA Lewis Research Center and the Air Force 
Wright Aeronautical Laboratories is also acknowledged. Finally the 
authors would like to thank Professor Toshi Kubota and Dr. Robert Gran for 
their suggestions regarding the interpretation of the experimental data. 
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ABSTRACT 


The investigatTo.i of an axi symmetric swirling confined jet flowfield, 
similar to the ones encountered in gas turbine combustors, has been car- 
ried out using a triple hot-wire probe. The raw data from the three sen- 
sors are digitized using ADC's and stored on a Tektronix 4051 computer. 
The data are further reduced on the computer to obtain time;-series for the 
three instantaneous velocity components in the flowfield. Finally, the 
time-mean velocities and the turbulence quantities are deduced. 

Qualification experiments were performed and where possible comparisons 
are made with the results of independent measurements. For example, the 
mean velocity components are compared with the results of five-hole pitot 
probe measurements. The major qualification experiments involved measure- 
ments performed in a non-swirling flow compared with conventional X-wire 
measurements. In the swirling flowfield advantages of the triple wire 
technique over the previously used multi-position single hot-wire method 
are noted. The results of the present measurements provide a data base 
with which the predictions of turbulence models in a recirculating swirl- 
ing flowfield can be evaluated in detail. 
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Calibration constants in Equation (3.1) 

Test section diameter 

Inlet nozzle diameter 

Hot-wire voltage 

Velocity normal to a single wire 

Velocity components in the coordinate system defined 
by a single wire. 

Velocity components in the coordinate system defined 
by the three wires. 

Velocity components in the coordinate system defined 
by the probe support axis. 

Instantaneous axial velocity in facility coordinate 
system. 

Instantaneous radial velocity in facility coordinate 
system. 

Instantaneous swirl velocity in facility coordinate 
system 

Pitcn factor. 

Yaw factor 

Number of points digitized per channel 

Radial location in the flowfield. 

Axial location in the flowfield. 

Effective cooling velocity acting on a wire. 

Side-wall expansion angle. 

Swirl vane angle. 

Probe pitch angle. 

Probe yaw angle. 


Refer to the three wires of a triple-wire probe. 
Root-mean-squared quantity. 


Time-mean averaged quantity. 
Fluctuating quantity. 
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1. INTRODUCTION 


1.1 Background and Previous Stadies 

Swirling confined jets have be^n under intensive investigation in the 
recent years. One application of such flowfields is a gas turbine com- 
bustor, an example of which is shown in Figure 1. Flowfields within such 
combustors typically have a rapid expansion and strong swirl imparted to 
the incoming air, which result in corner and central recirculation 
zones. The presence of swirl and the re-circulation zones increase the 
complexity in terms of fluid dynamic analyses of such flowfields. This 
complexity is further increased by the processes of combustion and heat 
transfer within the flowfield. Despite the complexity of combustor flows, 
significant progress is being made in their analysis as discussed in 
Reference 1. 

During an initial attempt to perform measurements of mean velocity and 
turbulence quantities, a much simplified combustor flowfield was model- 
ed. Figure 2 shows the schematic of idealized geometry with a control of 
side wall expansion and degree of swirl imparted to the flow. The figure 
shows various recirculation regions experienced in such a flowfield. This 
idealized swirling confined jet is being investigated at Oklahoma State 
University and at Dynamics Technology, with a variety of methods of 
approach. Analytically, a computer code (STARPIC) presents a theoretical 
model designed specifically to calculate the swirling confined jet flow- 
fields Experimentally, a series of flow visualization experiments 

coupled with 5-hole pitot probe measurements have been used to character- 
ize the time-mean flowfield^. 
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In addition to the mean properties, turbulence measurements are being per- 
formed using various hot-wire techniques. At Oklahoma State University, a 
technique of using a single, normal hot-wire at six different orientations 
has been developed which can measure time-mean velocities and normal and 
shear stresses. Some of the analysis of the measurements was performed at 
Dynamics Technology under the present sub-contract and is therefore 
included here as Appendix A. 

The six-orientation technique is basically a statistical technique applic- 
able to stati onary turbulent flows which uses mean voltage and root-mean- 
square voltage fluctuation measurements and an C'xtensive data reduction 
method to produce estimates of the components of mean velocity together 
with normal and shear turbulent stresses. 

1.2 Present Study 


At Dynamics Technology, a triple-wire measurement technique has been 
developed for the purpose of measuring mean velocity and turbulent 
stresses. The throe wires are simultaneously operated by three closely 
matched constant temperature anemometers, from which time series of these 
velocity components can be deduced. The instantaneous signals are digi- 
tized using ADC's which ctre controlled by a Tektronix 4051 computer. The 
subsequent data analysis is performed by digitally correlating the three 
signals and thereby obtaining the time-mean velocities and the velocity 
fluctuation products which can be averaged to produce Reynolds stress 
estimates. A significant advantage of the three-wire method is that the 
data reduction is performed on instantaneous measurements so that time 
records of velocity components (and fluctuations) can be obtained. Con- 
sequently more detailed information on the turbulence (including time 
domain information) can be obtained which is not possible with the single 
wire, six-position technique. 
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The data analysis is based on a study done by Yavuzkurt, CrawfoPu and Mof- 
fat ^ who used a triple-wire probe in a three dimensional turbulent flow- 
field. Yavuzkurt et al . developed a network of analog devices with which 
to process their signals. This data p.ocessing in the present study has 
been replaced by software coding operating on digitized signals. 

2. PROBE DESCRIPTION AND INSTRUMENTATION 


2.1 The Tri ple-Wi i g 


The three-wire probe used in the present study is DISA type 55F81. (See 
Figure 3.). The three sensors are 5 u-t’ diameter platinum plated tungsten 
wires. The wires are oriented in a mutually orthogonal array. The probe 
support lies symmetrically between the three wires making an angle of 54.7 
degrees with each of the wires or 35.3 degtvees with the plane defined by 
any of the two wires (see Figure 4). The gold-plated and stainless-steel 
prongs are embedded in three independent ceramic bodies mounted in the 
stainless-steel probe shank. 

The probe is operated by three separate but closely matched constant tem- 
perature anemometers (DISA type 55C12). The hot-wire raw signals are con- 
ditioned by DISA type 55N20 signal conditioners. The signal conditioners 
allow low pass filterin_ of the hot-wire signals at a desii'ed frequency. 

2 . 2 Experimental Fac ility 

Figure 5 shows the overall schematic of the entire facility. It consists 
of a 1400 C.F.M. fan manufactured by Joy Manufacturing Co. The fan drives 
the air througi a 9" diameter section of flow strai ghteners into a dif- 
fuser which has two perforated plates inside in order to prevent flow 
separation within the diffuser. The flow then passes through a turbulence 
management section which has several screens and flow straighteners to 
decrease the turbulence levels. The air then enters a contraction section 
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and exits through a round jet of 6 inch diameter and expands into the test 
section which consists of a five foot long plexiglass tube 12 inches in 
diameter. A swirler designed and built at Oklahoma State University can 
be installed at the nozzle exit along with the desired expansion block 
(see Figure 5.). 

2.3 The Probe Drive 


The probe drive, shown in Figure 6 was designed specifically for triple- 
wire measurements. The drive has three degrees of freedom: pitch and 

yaw, (necessary to aim the probe in the direction of the mean flow at any 
location in the flowfield) and vertical translation for surveying an axial 
station. 

2.4 Data Acquisit i on System 

The hot-wire signals from the three wires are digitized using analog-to 
digital converters which can scan each of the three channels simultane- 
ously up to a sampling rate of 9 kHz. The ADC's are controlled by the 
Tektronix 4051 laboratory computer. The digitized data are stored on mag- 
netic tape and further reduced and plotted using the laboratory com- 
puter. The entire data acquisition and reduction system is automated, 
thereby eliminating the intermediate handling of experimental data. 

3. THEORY AND DATA REDUCTTO* : 

3 . 1 Res:^»ise Equation s 

Before measurements in a turbulent flowfield, the triple-wire probe is 
calibrated in the uniform-flow region of the free jet. A calibration 
equation relating velocity to the hot-wire voltage for a single wire can 
be written as: 
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» A + B + C (3.1) 

where E is the hot-wire voltage and U,, is the component of velocity normal 
to the wire. For a triple-wire, three such calibrations are required for 
the three individual sensors. There are two methods of calibration 
reported in literature. One method, used by Gaulier ® for example, calls 
for orientation of the individual sensors normal to the calibration (in 
smooth potential flow) jet and thus obtaining an E vs curve fit given 
in Equation (3.1). The second method, used by Yavuzkurt et al . calls 
for orientation of the probe axis col inear with the cal i oration jet and 
resolving the velocity components in directions normal to the three wires 
to obtain an E vs 1)^ calibration for the three respective wires. For the 
present study, both the methods were tried and the method where the probe 
support axis is col inear with the jet was found to produce better results 
and hence is used in the analysis. Equation (3.1) can be now inverted and 
can be written in terms of the individual cooling velocities experienced 
by the three wires as: 


Uni = { -Bj + {Bj -4Cj(Aj - E^) 


Un^ = {-82 + {32 - 4C2(Aj - E^) 


(3.2) 


Un3 = (-83 + {83 - 4 C 3 (A 3 - E3) }^^^/2C3)^ 


Jorgensen^ introduced the sensitivity of a normal single hot-wire to flow- 
field angularity by an equation: 


z2 = 72 + Q2 72 + k2 w2 


(3.3) 
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where u,v, and w are velocity components In the coordinate system defined 
by the probe {see Flyure 7), z Is the apparent cooling velocity and the 
pitch and yaw sensitivity factors are defined by Jorgensen ^ as: 


G . JLiliijLJii 
u (v,w » 0) 

and 

K , V (u.»» * °) 

w (UoV = 0) 


(3.4) 


evaluated from calibrations perfonned in the three directions defined by 
Figure 7. The apparent cooling velocity z^ is the velocity that would be 
calculated from the hot-wire voltage (for the 1^^ wire) '“rom one of 

equations (3.2). 

As described in section 2, the wires of a triple-wire probe form a tri- 
hedral. Figure 8 shows the coordinate system (^q«>q»Zq) defined by the 
three wires. Equationf. (3.3) can now be written for each wire as: 


^2 + q2 w^ 

1 10 0 10 

z2 = G| v^ + w^ 

2 2 0 2 0 0 

2 ^2 ,2 
3 0 3 0 3 0 


(3.5) 


Equation (3.5) can also be written in the matrix form as: 
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G| K| 1 



1 
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.* ®3 '^3 J 


1 

0 K> 

1 


we can define a sensitivity matrix K given by; 


(3.6) 


K - 


Kf 1 

G| 

1 


'^3 


(3.7) 


Equation (3.5) can be Inverted to obtain velocity components In the Xg.yg 
and Zq directions (coincident with the orientation of the three wires); 




r ^ 



z? 

0 


1 


= K 

zi 

0 


2 

w2 


zi 

0 


L 3 J 


In calculating the velocity components Uq.Vq and Wq fron three squared 
values It is assumed that each component is the positive square root of 
the squared quantity. This is equivalent to assuming that the velocity 
vector lies within 54.7 degrees of the probe axis. In highly turbulent 
regions, points in time in which the velocity is outside the 54.7 degree 
cone are incorrectly resolved. This is perhaps the most serious source of 
error in the triple-wire measurement technique. 

In the calculation of the Uq.Vq and Wq velocity components at some j/Oints 
in time negative values of squared quantities are obtained (indicating 
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imaginary velocity components), te believe this is related to the situa- 
tion in which the direction of t.ie velocity vector is outside the 54.7 
degrees probe cone. There are two ways in which we handle such data. The 
first method is to completely discard the point from the data set. The 
second method is to set the individual velocity component (whose squared 
value is negative) to zero and include the data point in the set. Com- 
parison of final flowfield estimates calculated by the two ways provides a 
measure of the uncertainly of the quantity. 

Using the above described methods, the Uq, Vq and Wq velocity components 
are calculated at each point in time. These data can then be transformed 
into a coordinate frame x',y' and z' defined by the probe support as shown 
in Figure 4 by: 




COS45 COS35 

COS45 cos35 

COS45 COS35 


- ~ 
% 

V* 

= 

-COS45 

COS45 

0 


'^O 

w' 


cos45 sin35 

COS45 sin35 

-COS35 


/o. 


There is yet a final transformation that is perfonned to obtain the velo- 
city components In the facility coordinates (x,y,z). This transformation 
depends upon the orientation of probe axis in the flow. Assuming 0j , to 
be the pitch angle and 9^ be the yaw angle of the probe the final trans- 
formation equation can be written as: 


u 


cosSj^ cose^ -sin92 -cose^ sine^ 


u' 

w 

= 

cose^ sine^ cosf^ -sine^^ sine^ 


v' 

V 

• «■ 


sin6j 0 cose^ COS0J 

“ m 


w' 

L J 


(3.10) 
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where u Is the Instantaneous axial velocity, w instantaneous swirl 
velocity and v is the instantaneous radial velocity. Thus Equation (3.10) 
gives the time series for the three velocity components u,v, and w In the 
facility coordinate system x,y and z. 

Thus the time mean velocities are computed with 


n 



* I 
1>1 
N 

ii|/N 


. V 

1*1 

N 

v^/N 

(3.U) 

. V 

i, 

1«l 

Wj/N 



whe-'e N is the total number of points digitized per channel. 

Similarly, the turbulent normal stresses can be written as: 

N 

1 

U\ 



N 

(V, - v)‘ 

“y (3.12) 

N 

} { w . - w ) *' 

_ __ 
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Finally, the turbulent shear stresses can be written as 
N 

! - u)(v. - v) 


N 

I (u. - u}(w.- w) 

R 

N 

I (v j- v)(w. -w) 


Equation (3.2) through (3.13) are solved on a Tektronix 4051 computer 
using the digitized raw signals from the three wires. The computer is 
equipped with a graphics package which facilitates the plotting of radial 
and axial distributions of time-mean and turbulence quantities. These 
data are presented in the next section. 

4. RESULTS 


Measurements of time-mean and turbulence quantities using a triple-wire 
probe were carried out in the idealized, non-reactive flowfield at various 
axial locations. Experiments were performed with the swirler vane angle 
set at 38 degrees and the side wall expansion angle a of 90 degree. 

4.1 Data Sampling and Reduction ; 

As discussed earlier the signals from the three wires sometimes experience 
large differences generating negative values for squared quantities in 
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Equation 3.3. We suspect the problem to be caused by large local flow- 
field angles. Another possibility is that the turbulent scale sizes in 
the thin shear layer portions of the flowfield are small, and in some 
cases of comparable size to the triple hot-wire probe. In these cases 
probe resolution problems could lead to inaccuracies such as those 
experienced in the data. 

As discussed earlier data are either discarded or processed by setting 
individual unresolvable velocity components to zero. The effect of discard- 
ing data is difficult to estimate. Initially, a large number of data 
points are recorded so that enough points are retained to maintain statis- 
tical significance. This is tested by comparing the results of processing 
100, 200 and 400 data points and demonstrating "convergence" of the time- 
mean estimates. A second method involved comparing data calculated from 
the two methods of handling data with negative square roots discussed 
above. The most reliable method for estimating the accuracy of the data 
processing technique however, involved comparing results of the three-wire 
measurements with the results of independent measurements. In the swirl- 
ing flowfields comparisons are made with measurements of mean ''alocity 
components made with a five-hole pitot probe. In non-swirling flowfields 
comparisons are made with conventional x-wire measurements performed by 
Chaturvedi ^ of both mean and turbulent velocity components. Finally, 
comparisons are also made with turbulence measurements performed in an 
identical swirling confined jet flowfield ^ using a six-orientation, 
single hot-wire technique. However, it appears that the present three- 
wire probe technique is in fact more accurate ana reliable than the six- 
orientation, single wire method. This will be discussed further with the 
presentation of the measurement results. 
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4.2 Techniqufc Verification in Non-Swirling Floi ; 

Preliminary measurements were performed in a non-swirling confined jet at 
x/D = 2.0 with (^ = 0 and a = 90 degrees. The results are compared with a 
similar study performed by Chaturvedi^ in a corresponding flow situation 
using a cross-wire probe. Figure 9 shows the radial distribution of time- 
mean and turbulence quantities. The two studies show very favorable 
agreement. Overall, the results show impro''^.ment over the results obtain- 
ed by the six-orientation hot-wire technique reported in Reference 4, 
Appendix A. In particular, the radial turbulence intensity, when measured 
by the six-orientation technique was lower than the one measured by Cha- 
turvedi^ (at the x/0=2 station; see Reference 4). The triple-wire 
measurements show the radial turbulence intensity to agree very well with 
the Chaturvedi measurements. Generally speaking, the non-swirling con- 
fined jet measurements verify the triple-wire technique. 

4.3 Mean-Velocity Field in Swirling Flow ; 

In a swirling confined jet flowfield, the axial and azimuthal velocity 
components are the two dominating components. The swirling confined jet 
flowfield measurements reported here have all been performed with para- 
meters D/d=2, the swirl angle $=38 deg and the expansion angle Q!=90 deg. 
These parameters all correspond exactly to those used by Yoon^ in an 
experimental study of the mean flow of the confined swirling jet. Yoon's 
measurements were all performed with a five-hole pitot probe at axial 
coordinates x/D=0.5,1.0,1.5, etc. In addition, he reported radial distri- 
butions of mean axial and swirl velocity components in the confined jet 
exit plane (at x/D=0). These velocity measurements indicate that the cen- 
tral recirculation zone starts at the base of the center hub. Outside of 
the center hub wake, the angle of the flow is approximately 28 deg, to the 
facility axis (compared with the blade angle ^ of 38 deg). Figure 10, 11 
and 12 show the radial profiles of the axial, radial, and azimuthal com- 
ponents of mean velocity. The data are compared with the five-hole pi toe 
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probe measurements^. The two studies are in excellent agreement in the 
case of the axial component of mean velocity. The radial component of 
mean velocity also shows favorable agreement with the five hole pitot 
probe measurements for x/D 21. The swirl component measured by the 
triple-wire however, is up to twenty five percent lowe. than the one 
measured by the five- hole pitot probe at axial location of 1.5 and 
higher. At lower axial locations, the two measurements seem to agree very 
well . 

As mentioned earlier, the measurements have been performed in a swirling 
confined jet in which the angles of the swirl blades have been set to 
38°. After expanding to the larger diameter (D=30 cm) from the exit dia- 
meter (D=15 cm), local flow angles significantly increase to near 60° to 
the facility axis. The differences in the swirl velocity component 
measurements, seen in Figure 12, represent differences in local flowfield 
angles of up to 10°. The reason for this discrepancy is not known at this 
time. It is important to note however, that measurements of the swirl 
velocity component made with both the present three-wire method, and with 
the six-orientation single-wire technique are in almost perfect agreement 
(Compare solid symbols in Figure 12 with solid symbols in Figure 10, 
Appendix A). 

4.4 Turbulence Intensity Measurements : 

The primary purpose of the hot-wire probe is to measure the velocity fluc- 
tuations in a flowfield, taking advantage of the inherent fast time 
response of the instrument. A triple-wire probe makes it possible to 
measure the three components of the turbulence velocity in the swirling 
flow under investigation. Figures 13, 14 and 15 shows the radial distri- 
butions of the mean-square of three components of the turbulence velo- 
city. The triple-wire data are compared with data obtained with the six- 
orientation hot-wire technique^. The two studies show reasonable agree- 
ment in quantifying the mean square of three components of the turbulence 
vel oci ty. 
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In general the turbulent fluctuations are higher in regions of higher mean 
shear, such as evidenced on the edge of the recirculation region at an 
axial station x/0=0.5. However, there appear to be substantial turbulence 
levels throughout the recirculation region which we suspect is caused by 
the large scale "breathing" of the region. Understanding of the phenome- 
non will require more detailed measurements and signal analysis than was 
possible in this study. 

4.5 Turbulent Shear Stresses ; 

Shear stresses are the most difficult turbulent quantities to measure 
accurately in a complex three-dimensional flowfield. Uncertainities in 
measurements of mean velocities and turbulence intensities are multiplied 
in determination of the shear stresses. Figures 16, 17 and 18 show the 
turbulent shear stresses deduced from the three-wire hot-wire measure- 
ments. When compared with the corresoonding measurements performed in the 
non-swirling confined jet, (Figure 9), it can be seen that at most axial 
stations the radial location of highest turbulent shear stresses has moved 
outwards, corresponding to the shear layer displacement induced by the 
centrirugal forces. Note, that the measurements produce primarily 
negative values of the shear stress component u'v' which is expected in 
the viscinity of the recirculation region where the gradient of the axial 
velocity (5^/ar) is predominantly positive. 

4.6 Comparison of the Three-Wire and the Six-Orientation Single 
Hot-Wire 1-lea sureinent Techniques : 

The three-wire measurement technique was adopted to the present research 
study with the idea that improved experimental data would result. In par- 
ticular, accurate estimates of turbulent shear stresses are difficult to 
obtain with any measurement technique (laser Doppler and hot-wire anemo- 
metry being the two most practical). Measurements were previously per- 
formed in the Oklahoma State University confined jet facility (with 
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a-90 deg and $=38 deg) with the six-orientation single hot-wire tech- 
nique. The resulting data are shown on many of the present figures and in 
more complete form in Reference 4 included as Appendix A. 

The u'v' data deduced from the six-orientation single hot-wire measure- 
ments are all shown as positive because the method used in that study was 
unable to distinguish between positive and negative radial velocity com- 
ponents, If the six-orientation single-wire data were assigned their 
expected sign {to correspond to the results of the three-wire measure- 
ments) then the two sets of data would be in closer agreement. However, 
the u'v' data from the six-orientation single wire technique are usually 
larger in magnitude than the triple-wire data and exhibit significantly 
more scatter. 


It could be argued that the u'v' data from the present method are more 
accurate than the single wire measurements. First, the time mean radial 
and axial velocity components are more accurately determined (as shown by 
the comparison with the 5 hole pitot probe measurements). Second, the 
turbulence measurements made in the non-swirling flow are slightly more 
accurate (as shown by comparison with Chaturvedi's x-wire measure- 
ments^). Finally, the magnitudes of all the turbulent stresses, as deter- 
mined by the three-wire method, are of similar magnitude to the corres- 
ponding data measured in the non-swirling flow, if account is taken of the 
large changes in the recirculation geometry. The six-orientation single 
wire measurements indicate a larger u'v' component of the turbulent stress 
in comparison with the non-swirling confined jet. This does not seem 
reasonable to us because we would expect the centrifugal field to stabi- 
lize the turbulence rather than amplify it. Thus, in summary, the present 
three-wire measurements appear to be producing more reliable results than 
the previous single wi re-measurements. 
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Uncertainty estimates were made on the data by two methods. First the 
sensitivity of the calculations to small changes In primary measured quan- 
tities was determined. These calculations showed significantly less error 
magnification than was determined for the six-orientation single hot-wire 
method^. Additionally, hot-wire signals were analyzed using the two met- 
hods of handling data with negative squared quantities discussed earlier 
In Section 3.1. The difference between data determined In these two man- 
ners provides a measure of the uncertanity In the final estimates. A 
representative sample of the data Is shown In Figure 19. Included are the 
profiles of mean axial velocity, the axial component of turbulence velo- 
city and the u' v* component of the turbulent shear stress, measured at 
x/D=1.0. The solid circles represent estimates in which negative squared 
quantity data are eliminated from the data set. The open circles are 
estimates in which individual velocity components (in wire coordinates) 
are set equal to zero when encountering a negative squared quantity. 
These two estimates define an anticipated bound for the actual data. 

Included on the figures (with the x symbol) are the estimates obtained 
from the slx-orientatiun single hot-wire measurements^. Note that there 
Is very good agreement between the two methods in the estimates of axial 
mean velocity and axial turbulence velocity. For this comparison, u'v' 
values measured by the six-orientation technique have been assigned the 
appropriate sign corresponding to the results of the triple-wire study. 
In the turbulent shear stress estimates, there is typically more scatter 
In the single wire measurements with indicated magnitudes usually larger 
than the triple-wire measurements. 

A final point to be made regarding the comparison between the two hot-wire 
methods concerns the ability of the triple-wire technique to produce time 
series data as opposed to the single-wire technique whose reduced data is 
always ’n time averaged form. A detailed discussion on time-series data 
is presented in the following section. 
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4.7 Real-time Measurements ; 

The triple-wire study presented here is unique in the sense that it pro- 
vides time-series for the velocity vector sensed by the three wire 
probe. The previously used multi-orientation single hot-wire technique 
uses the averaged and root-mean-squared fluctuation values of hot-wire 
voltage for data reduction. The triple-wire data acquisition is automated 
and signals are digitized using rer.l-time analog to digital converters. 
The triple-wire technique thus has a definite advantage over the six- 
orientation single wire method. 

Figure 20 shows the time-series of the velocity experienced by one of the 
three wires. This record was obtained in the swirling confined jet at 
x/D=1.0, r/D=0.3 with (t»=38 degrees, and a=90 degrees. The figure shows 
large variations in velocity occuring at frequencies up to about 200 Hz. 
Most of the turbulence energy is carried at these low frequencies. In 
addition, there are large peaks occuring at even lower frequencies such as 
20 Hz. 

Figure 21 shows the power spectral density of the signal from a single 
hot-wire in the confined jet flowfield. In this record 1024 points were 
digitized at a rate of 32 kHz. The plot shows that most of the energy of 
turbulence is contained in frequencies lower than 1 kHz. A decision was 
then made to set tiie data digitization rate at 2.0 kHz for production run 
experiments. 

A review of Figure 20 and 21 suggests that the hot-wire is sensing a con- 
siderable amount of low frequency, undoubtly large scale fluid motions. 
It is possible that these large scale motions are coupling with the acous- 
tic modes of the test chamber. For example, the fundamental plane wave 
mode for the test chamber is approximately 45 Hz and the spinning modes 
are correspondingly higher in frequency. These acoustic modes surely 
cause the recirculation regions to "breathe" producing considerable energy 
in low frequency fluctuations. This low frequency energy is mostly 
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deterministic in nature and therefore is not best represented by 
conventional turbulence models. The triple-wire hot-wire technique 
developed in this study could be used to conduct a detailed study of the 
large scale motions and their acoustic interactions in the swirling 
confined jet flowfield. Such a study would provide valuable information 
to aid in the interpretation of combustion instabilities in axisyn«netric 
combust. '^rs. 
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5. CONCLUSIONS 


5.1 Summary 

The measurements of time-mean velocities and shear stresses in a complex 
three-dimensional flowfield have been carried out using a triple-wire 
probe. The salient features of the outcome of the present study can be 
buiMMiari zed as follows: 

1. Mean velocities measured by the triple-wire probe are in very good 
agreement with measurements perfonned with a five-hole pitot probe in 
an identical flowfield, 

2. Ii; the non-swirling confined jet, turbulence measurements are in 
very good agreement with previous measurements performed by Chatur- 
vedi^ using a conventional x-wire probe. 

3. In the swirling confined jet flowfield the turbulence intensities 
and turbulent shear stress estimates ootained from the triple-wire 
technique are in general agreement with the corresponding measurements 
performed with the previous six-orientation single hot-wire method. 
In fact, the present measurements appear to be more accurate (and 
hence more reliable) insofar as they retain the time series 
information lost in the single-wire method. 

5.2 Further Study : 

In the present study the triple-wire measurements have been performed for 
a confined jet flowfield with the swirl vane angle ^38 deg. with an 
expansion angle a=90 deg. The measurements have produced detailed baseline 
turbulence data which can be used to evaluate the applicability of various 
turbulence models in this swirMng flowfield. 
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An obvious extension of this work would Include a parametric study of the 
effect of vane angle and expansion angle on the turbulence properties. 
The ability of the triple-wire probe to perform real-time measurements In 
complex (swirling) flowflelds suggests a detailed examination of velocity 
transients and their coupling with acoustic modes. This type of study 
should also be performed with downstream blockage found In realistic com- 
bustor geometries. 
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Figure 1. Typical Axisymmetric Combustion Chamber of a Gas Turbine Engine 
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CRZ - Corner Recirculation Zone 

CTRZ - Central Toroidal Recirculation Zone 


Figure 2. Idealized Combustor Flowfield. 
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Figure 5. Schematic of Overall Facility. 
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Figure 6. The Three Dimensional Probe Drive. 
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Confined Jet with #=38 and ^=90 
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Figure 18. Radial Distribution of Shear Stress v'w' in Swirling Confined Jet 
with *=38° and a=90°. 
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Figure 19. Bands of Uncertainities in Measurements of Mean and Turbulence Quantities 
Measured at x/D = 1.0 with *=38° and 90°. 
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Figure 20. Time-Series of Velocity Experienced by one of the Wires in Swirling Flow 
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Figure 21. Power Spectral Density Plot for Non-Swirling Confined Jet Flow, 
at x/D = 1.0 and r/D = 0, with <t>= 0°. 
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APPENDIX A 

Turbulence Measurement In A Confined Jet 
Using A Six-Orientation Hot-Wire Probe Technique 
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Abstract 


The s1 x-orlentatlon single hot-wire technique 
has been applied to the complex flowfleld of a 
swirling, confined jet. This flowfleld, which 
contains a rapid expansion with resulting recircu- 
lation regions. Is typical of those found In gas 
turbine engines and ramjet combustors. The pre- 
sent study focusses on turbulence measurements In 
such a flowfleld In the absence of chemical 
reaction. 

A modification to the si x-orlentatlon hot-wire 
technique developed by King has been made, which 
incorporates the determination of turbulent shear 
stresses (in addition to normal stresses) and 
ensemble averaging of redundant turbulence output 
quantities. With this technique, flowflel d surveys 
have been performed in both swirling and nonswlrl- 
ing axl symmetric confined jets. Where Independent 
data exist, comparisons have been made which 
demonstrate the reliability of the technique. 
Finally, a sensitivity analysis of the data reduc- 
tion technique has been completed which forms the 
major Ingredient In an uncertainty analysis. 

Nomenclature 


A,8,C Calibration constants in Equation 1 
AO, BO, CO Cooling velocity functions In Table 1 
0 Test section diameter 

d Inlet nozzle diameter 

E Hot-wire voltage 

II Velocity function for axial velocity 

W Velocity function for azimuthal velocity 

V Velocity function for radial velocity 

G Pitch factor 

K Yaw factor 

K, , Covariance for cooling velocities 

Zp, and Zq 

P,0,R Selected hot-wire probe positions 

Re^ Inlet Reynolds number 

u .Axial velocity I in 

V Radial velocity > coordinates 

w Azimuthal (swirl) velocity ) on facility 

u,w,v Three components of velocity In probe 

coordinates defined by Figure 5 
x.r.e Axial, radial, azimuthal cylindrical 

polar coordinates 

Z Effective cooling velocity acting on a 

w1 re 
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** Senior Research Scientist, Member AIAA 
t Graduate Student, Student Member AIAA 
tt Professor, Associate Fel ow AIAA 

( <ip>ri|thi ~ Amriiraii Ibililul. of Arraniulirt aiHi 
A.lriinaultc.. lac.. 19S2. All riihu mrr.rd. 


a Side-wall expansion angle 

I Correlative’’ coefficient (estimated) 

1 j between cooling velocities of 

adjacent wire orientations 
Variance of a given o itity 
♦ Inverse function of calibration equation 

i Swirl vane angle 

Subscripts 
1.2. 1,4, 5, 6 

1.j 
P.Q.R 
rms 

Superscripts 

Tlme*mean average 
' Fluctuating quantity 

1. Introduction 

1.1 The Gas Turbine Combustor Flowfleld 


Refers to the six probe measuring 
positions 

Dumihy Indicles which take the values 
1 to 3 

Refers to the three selected 
measuring positions 
Root-mean- squared quantity 


Recent emphasis on fuel economy and pollutant 
suppression has sparked a renewed Interest in gas 
turbine combustor analysis. A typical axl sym- 
metric gas turbine engine combustor Is shown In 
Figure 1. Flowflelds within such combustors typi- 
cally have a rapid expansion and strong swirl 1m- 
,/drted to the Incoming air, which result In corner 
and central recirculation regions. The swirling, 
recirculating, turbulent flows within combustors 
present one of the more difficult fluid dynamic 
problems to analyze. This complexity Is Increased 
many fold by the processes of combustion and heat 
transfer within the flowfleld. Despite the com- 
plexity of combustor flows, significant progress 
is being made In their analysis.^ 


Viin«) S*Kl*r 



Figure 1. AkI symmetric Combustor 
of a Gas Turbine Engine 
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Tne present paper reports on research which is 
part of an extensive experimental and computa- 
tional study of gas turbine flowfields In the 
absence of combustion. Figure 2 shows the charac- 
teristics of the simplified flowfield being inves- 
tigated. Flow enters through a Jet of diameter d 
into a tube of diameter 0, after being expanded 
through an angle o. Before entering the tube, the 
flow may be swirled by a swirl er located upstream 
of the inlet plane. Shawn schematically are the 
corner recirculation zone (CRZ) and the central 
toroidal recirculation zone (CTRZ) which are typi- 
cally present In these flows. 



Figure 2. Idealized Combustor Flowfield 


The swirling confined jet flowfield shown i > 
Figure 1 is being Investigated at Oklahoma State 
University and at Dynamics Technology, with vari- 
ous methods of approach. Analytically, a computer 
program (STARPIC) has been assembled which is de- 
signed specifically to calculate the swirling con- 
fined jet flowfields.^ Experimentally, a series 
of flow visualization experiments coupled with 5- 
hole pitot probe measurements have been used to 
characterize the time-mean flowfield.^ Hot-wire 
measurements of the turbulence properties are also 
being conducted. This paper reports on the Ini- 
tial results of the hot-wire measurements in the 
confined jet flowfield. 

Several studies on time-mean flowfields of the 
type just described have been carried out usings 
various turbulence measuring techniques.'*"^ ^ 
Unfortunately, most of the techniques used do not 
give complcce and detailed Information about the 
flow in terms of all Its time-mean and turbulence 
quantities. In addition, no experiments have been 
performed on the specific geometry of the present 
study in the presence of inlet swirl. To develop 
further the flowfield computational techniques, 
including the turbulence modeling, there is a 
strong need to obtain experimental estimates of 
the tii'^bulence and mean flow quantities in such 
flows. 

1 .2 The Turbulence Measurement Problem 

Turbulence mea .urement in a complex flowfield 
has always been a complicated problem encountered 
by engineers. In the past, turbulence phenomena 
have beei’ Mscussed by various authors in detail 
and varic,^ methods of turbulence measurements 


have been suggested. One of the most widely 
used instruments to obtain turbulence quantities 
is the hot-wire anemometer, the most common of 
which is a single hot-wire. When used at a single 
orientation and in a two-dimensional flow with a 
predominant flow direction, a single hot-wire can 
measure the streamwise components of the time-mean 
velocity and the root-mean- square velocity fluc- 
tuation at a particular location in the flowfield. 
A two-wire probe can be used to determine the 
time-mean velocities, streamwise and cross stream 
turbulence Intensities, and the cross correlation 
between the two components of the velocity fluc- 
tuations.^**'!^ 

Hot-wire measurements in complex three-dimen- 
sional flowfields are considerably more difficult 
than in one- or two-dimensional flowfields in 
which the mean flow is predominantly in one direc- 
tion. To measure the three velocities and their 
corresponding fluctuating components in a three- 
dimensional flowfield such as encountered in com- 
bustor simulators, there are two methods tnat can 
be employed at a point in the flowfield: 

1) A three-wire probe used with a 
single orientation. 

2) A single- or double-wire probe 
used with multi-orientation. 

The three-wire probe technique permits the 
necessary simultaneous measurements from which 
three Instantaneous velocity components can be 
determined. The appropriate signal processing can 
produce estimates of mean velocity components and 
normal and shear turbulent stresses {such as 

u' ’ and u‘ v' 1 . 

The three-hot-wire probe technique is signifi- 
cantly more complex than the single wire multi - 
orientation techniques. A multi-dimensional probe 
drive is required to orient the probe in epproxi- 
mately the mean flow direction. Also, sophistica- 
ted signal processing electronics is required to 
handle the three instantaneous hot-wire voltages. 
Finally, the three-wire probe typically has less 
spatial resolution in comparison with a single 
w1 re probe. 

Multi-orientation of a single hot-wire is a 
novel way to measure the three components of a 
velocity vector and their fluctuating components. 
A method devised by Dvorak and Syred^* uses a 
single normal hot-wire oriented at three different 
positions such that the center one is separated by 
45 degrees from the other two. The velocity 
vector at a 1 ation is related to the three 
orthogonal components using pitch and yaw factors 
as defined by Jorgensen. The data are obtained 
in the form of mean and root-cean-square voltages 
at each orientation. However, the measurements 
done with a single wire do not supply all the 
information needed to obtain the turbulence quan- 
tities. Therefore, in addition to a single wire, 
Dvorak and Syred used a cross-wire probe to obtain 
the covariances between the voltages obtained at 
adjacent hol-wire orientations. 

King^® modified the technique developed by 
Dvorak and Syred. Hi? method cells for a normal 
hot-wire to be oriented th ough six different po- 
sitions, each orientation separated by 30 degrees 
from the adjacent one. Mean and root-mean-square 




voltages are measured at each orientation. The 
data reduction Is performed using some assumptions 
regarding the statistical nature of turbulence, 
making It possible to solve for the three time- 
mean velocities, the three normal turbulent 
stresses, and the three turbulent shear stresses. 

1.3 The Scope of the Present Study 

In the present study, the si x-orlentatlon 
single normal hot-wire technique Is being employed 
to obtain the turbulence quantities In the combus- 
tor simulator confined jet flowfleld. Measure- 
ments have been carried out for both swirling and 
nonswlrllng flow with axpanslon angles of 90 
degrees (sudden expansion) and 45 degrees (gradual 
expansion). Only the 90 degree angle data are 
presented here and the Reynolds number of the 
Inlet flow 1s * 5 x 10'' which Is comparable with 
aircraft combustor flows (although our experiments 
are performed In nonreacting flows). The data 
reduction procedure extends King's technique to 
obtain turbulent shear and normal stresses using 
six basic response equations representing the six 
orientations of a normal hot-wire positioned in 
the flowfleld. Certain modifications are made In 
the procedure to calculate covariances which are 
an lnteg-al part of the data reduction procedure. 
An uncertainty analysis Is performed on the tech- 
nique which reveals the sensitivity of this tech- 
nique to various Input parameters discussed In the 
later parts of this paper. Some of the turbulence 
quantities obtained are compared with measurements 
performed by Chaturvedl^ using a crossed-wire 
probe In a corresponding flow situation. 


2 . Experimental Facility and Instrumentation 
2.1 Idealized Flowfleld 

The facility, designed and built at Oklahoma 
State University, Is a simulation of a typical 
a x1 symmetric combustion chamber of a gas turbine 
engine shown In Figure 1. The schematic of the 
test facility with the Idealized flowfleld Is 
shown In Figure 3. Anbieit air enters the low- 
speed wind tunnel through a foam air filter. The 
air then flows through an axial flow fan dv'lven by 
a 5 h.p. varl drive motor. Thus, the flow rate can 
be varied for different test conditions. The flow 
passes through a turbulence management section 
which has two f1ne-mesh screens, a 12.7 cm length 
of packed straws, and five more f1ne-mesh screens. 



Figur* 3. Schematic of the Experimental 
Facility 


The axl symmetric nozzle was designed to pro- 
duce a minimum adverse pressure gradient on ' le 
boundary layer to avoid flow unsteadiness asso- 
ciated with local separation regions. The area 
ratio of the cross sections of the turbulence 
management section to that of th; nozzle throat Is 
approximately 22.5. The diameter, d, of the 

nozzle throat 1$ approximately lo cm. 

The test section is composed of a swirl er 
(optional), an expansion block, and a long plexi- 
glass tube. The expansion block, attached after 
the swirler. Is a 30 cm diameter annular disk of 
wood. At present, there are three expansion 
blocks, and the appropriate choice gives a « 90, 
70, or 45 degrees. The flow Is expanded Into a 
plexiglass tube of diameter, 0, of 30 cm. thus 
giving a diameter expansion ratio (D/d) of 2. The 
test chamber has no film cooling holes or dilution 
air holes, and the chamber wall of the test sec- 
tion 1s a constant diameter pipe. 

2.2 Hot-Wire Instrumentation 

Jhe anemometer used for the present study 1 s 
OISA type 55M01, CTA standard bridge. A normal 
hot-wire proba. OISA type 55P01, Is used In the 
experiments. This probe has two prongs set 
approximately 3 im apart which support a 5 m<’< 
diameter wire which 1s gold plated near the prongs 
to reduce end effects and strengthen the wire. 
The mean voltage Is measured with a Hickok Digital 
Systems. Model DPIOO. Integrating voltmeter and 
the root-mean-square voltage fluctuation Is mea- 
sured using a l^wlett Packard. Model 400 Ilk. AC 
voltmeter. 

The hot-wire Is supported In the facility b> a 
traversing mechanism shown schematically In Figure 
4. It consists of a base that Is modified to 
mount on the plexiglass tube of the test section 
at various axial locations. T^e hot-wire probe Is 
Inserted Into the tube through a rotary vernier 
and the base. The rotary vernier Is attached to a 
slide which can traverse across the flow chamber. 
Thus. It Is possible for the probe to be traversed 
to any radial location at selected downstream 
locations In the flowfleld and to be rotated 
through 180 degrees. 
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Figure 4. Mot-wire Probe Mounted on the 
Test Section 

* Provided for Information and not necessarily a 
product endorsement. 
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2.3 Calibration Nozzle 
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The hot-wlrr is calibrated In a small air Jet. 
The facility consists of a compressed air line, 
which delivers the desired flow rate through a 
small pressure regulator and a Fisher and Porter 
Model IOA1735A rotameter. The jet housing con- 
sists of an effective flow management section fol- 
lowed by a contoured nozzle with a 3.5 cm diameter 
throat. A rotary table is used to hold the probe 
while It 1$ being calibrated In three different 
orientations. 

3. Hot-W ir e Data Analysi» 

3.1 Hot-Wire Response Equations 

The six-orientation hot-wire technique 
requires a single, straight, hot-wire to be cali- 
brated for three different flow directions 1n 
order to determine the directional sensitivity of 
such J probe. The three directions and three 

typical calibration curves are shown In Figure 5. 
In these relations, tildes signify components of 
the Instantaneous velocity vector In coordinates 
on the probe. Each of the three calibration 

curves Is obtained with zero velocity In the other 
two directions. The calibration curves demon- 
trate that the hot-wire Is most efficiently 
cooled when the flow Is in the direction of the u 
component, whereas, the wire Is most Inefficiently 
cooled when the flow Is In the direction of 
the w component. Each of the calibration curves 
follows a second order, least square fit of the 
form: 




A^ . 




n) 


where G and K are the pitch and yaw factors 
defined by Jorgensen^ ^ 

- v(u and w»0) 
u ■ * Z ' " ' 

u(y and w*0) 

K « v(u and w«0) 
w(u and v*0) 

which are evaluated from the three calibration 
curves (Figure 5) for a constant value cf E^. 
Equation (3) shows that ^he pitch and yaw factors 
a^e calculated with the v component 1 ■ 2 in equa- 
tion (1) of the effectiv cooling velocity as the 
relrrence. Ther^i' . iC calibration constant'.'' 
usedjn equation (i; •: -e the coefficients In the E 
vs. v calibration of Figure 5., I.e., In a general 
flowfleld: 


WV VW • 


, and 


(3} 


E2 X A^ + ♦ C.^ Z 

with Z as given In Equation (2) above. 

Figure 6 shows the pitch and yaw factors as a 
function of hot-wire voltage determined from the 
calibration curve of Figure E. Both factors v».ry 
with hot-wire voltage, but the yaw factor Is far 
more sensitive. The si;nslt1vlty analysis dis- 
cussed In the next se;t1on demonstrates that 
uncertainties associated with the varying pitch 
and yaw factors do not sirlously affect the accur- 
acy of the estimated flo«j quantities. 


which Is an extinslon of the commonly used Kind's 
law. In this equation, ^A^ , Bj, and are 
calibration constants and u. tan take on a value 
of u, V, and w for the three' calibration curves, 
respecti vely. 



Figures. 1 he Three-01 recti onal Mot-Wlre 
Cellbratlon 


When the wire Is placed In three dimensional 
flowfleld, the effective cooling veinc! y experi- 
enced by the hot-wire Is: 



iC } 1 12 3 i 14 M 

NOT -Witt voLiAOi fivoUi; 


Figured. Pitch and Yaw Factors Plotted Against 
Hot-Wire Mean Effective Voltage 


■•■o carry out measurements In the confined jet 
flowfleld, the wire Is aligned In the flow In such 
1 vay that In the first orientation, the wire Is 
-ipmul to the flow In the axial direction and the 
prob< voordl nates coincide with the coordinates of 
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the experimental facility. Thus, the six equa- 
tions for the instantaneous cooling velocities at 
the six orientations, as given by Klng^® are: 


probe. Therefore, Equations (5) must be expressed 
In terms of mean and root-mean-square values. 
Equation (1) can be written ns; 


= v^ + (j^u^ + 

Z^ v^ + G‘(u cos 30®+ w sin 30®)^ 
+ k2(w cos 30®- u sin 30®)^ 

Z| = v^ + G'^lu cos 60®+ w sin 60®)^ 
+ K^(w cos 60®- u sin 60®)^ 


Z,= [-B2+{b2+ 4C2(«2-E2)}^''^/2C2]^ . (6) 

The above equation Is In terms of Instantaneous 
velocity Z^ and 1 nst; itaneous voltage E^ . In 
order to obtain an expression for time-mean velo- 
city as a function of time-mean voltage, a Taylor 
series expansion of Equation (6) can be carried 
out as follows: 


= v^ + G^w^ + k^u^ (4) 

Z2 = v^ + g 2 (w sin 120®+ u cos 120®)2 
+ k 2 (u sin 120®- w cos 120®)^ 

Z| = v2 + g2(w s1n 150®+ u cos 150®)2 
+ K^u sin 150°- w cos 150°)2 


Zi= Z^(E^+EJ)= (|»[E^)+ ^||_+ 




17) 


where ^ 


The Taylor series Is truncated after second order 
terms assuming the higher order terms to be rela- 
tively small. Time averaging both sides of the 
above equation and employing the fact that E'=0, 
yields : 


Solving simultaneously any three adjacent equa- 
tions provide expressions for the instantaneous 
values of the three velocity components, u, w, and 
V in terms of the equivalent cooling velocities 
{Z^, Z2 and Zj for example, when the first three 
equations are chosen). Thus, the general form of 
the Instantaneous velocity components Is given as: 


U = [(AO + [AQ2 



1 

(G2-K2) 


1/2 

] 


« n 2 1/2 . 1/2 

W = [ ( -AO + (AQ2 + } i ] (5) 

^ (G2-K2) 


V = [ CO - (A02 + ^ 

(G2-K2) ^ 

The values of AO, BO and CO depend on the set 
of the three equations chosen and are given 1n 
Table 1 for appropriate equat on sets. 


Zj - 4 + 

^ ^ SE2 


( 8 ) 


To obtain 1'.'^ - a?. , the relationship as given by 

Hlnze'^ Is:’ ^ 

Z^^ ■ - (Expec [Z.])2 . (9) 

Using Equation (8) as the basis, Expec [Z^ 
and (Expec[Z^3)2 can be evaluated and substituted 
Into Equation (8) to get: 


z;2 . 




2 f 1 92^ 

9rr®E,- ^ ® 


■1 “I 




( 10 ) 


Thus, Equations (8) and (10) give the mean and 
variance of effective cooling velocities In terms 
of the mean and variance of the appropriate vol- 
tages. 


IARLE 1 


values r>f AO, M, anrt CO in Vtrioot Equation 


'Ijuatlon set 
P. '.R Cfioice 

1 

AO 

BO 

CO 

1. 1 

- //I 

. J 

.-2Z‘ ♦ 3Z2 - Z2) 
; 2 3 

[li - 4 Z2) 

' l 2 3 

Z. 5. 1 

' 7*' - ' 

:-l‘ . 17- - 7l>] 

(Z7 - Z> . 77 ; 


; i 

i i •> 

2 J % 

J. 4. 

; . 71‘ , Zi 1 

!77 - zn 

(Z7 - V . Z7) 


1 «• 5 

) 5 

i a i' 

*. s. ft 

■l- • z-'i 

1-777 , . .7\ 

(77 - 77 ♦ 7») 


S a 

1. S 6 

a 5 ♦ 

S, k, \ 

» V' 

1 - 7 ’ • Jz' - 777 ' 

(77 - Z7 . 7') 

» 4 ) 


'1 a 

5 6 : 

ft. 1. ^ 

:-7' . 21‘ - zn 

(-17 . Z7] 

(Z7 - 77 . 7») 


« 1 2 

« 2 

« 1 2 


In a 3-d1mens1onal flow. It Is usually desir- 
able to obtain the mean and variance for the Indi- 
vidual velocity components In axial, azimuthal, 
and radial directions, and also their cross cor- 
relations. The procedure to obtain the mean and 
variance of the Individual velocity components Is 
the same as for the effective cooling velocities 
except that u, w and v are functions of three 
random variables and there are extra terms In the 
Taylor expansion to account for the covariances of 
the cooling velocities. Thus, the axial mean 
velocity component as given by Dvorak and Syred,^* 
and K1 ngi® Is: 


=U(Zp,Zo,Z„] 


1 y 

1*1 


d2u 

9Z2 


3 

KJ 


92y 

'SI. 


1^ 


z z 


( 11 ) 


However, these equations cannot be directly k „ . 

used because 1t Is Impossible to obtain Z^, Z^ and +'nere l^z^Zj covariance of the cooling velo- 

Z3 at a single Instant 1n time with a single wire cities Z^ and Zj and Is defined as: 


5 
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1 T - 

z * - ^ <it (12) 

Identical expressions for w and v can also be 
obtained In terms of W and V, respectively. 
Derivatives of the form a^U/aZ^.bZj are determined 

analytically from equations (5) and Table 1. 


tions would be such that their contribution to the 
cooling of the wire would be related by the cosine 
of the angle between the wires. This assumption 
leads to the following three valuos of the corre- 
lation coefficients. 



= cos 30 = 


0.9 


(16) 


Also, the normal stresses are given as: 


3 

11*2 = y 


■ dU 


2 . V r aU ou „ 

5 ST • »Z 7 • "z,2 


au 


1=1 ""1 






Y7 7 = cos 30 


0,9 


To relate Y7 7 with y? 7 andy7 7 , 

'ZpZR 'ZpZq 

duced the following relationship: 


King 1ntro- 


1 V _2 . V 

f 1 — - I T T "S T ~ * ^Z Z 

^ 1=1 az| ^1 Kj 


(13) 


with similar expressions for w'^ and v'^ . 

Finally, the expressions for shear stresses as 
given by Dvorak ard Syred^*-' are of the form: 

-tt au bV 2 ^ ? V au av ^ 


r 1 V . V r a 2 u 


1=1 az 2 1 i^^jj 1 j 1 j 


%Zp = ^%2q%Zr 

where n Is given a value of 0.8, Hence z 
becomes: ^ ^ 

= (0.8)(0.9)(0.9) = 0.65 (18) 


The three covariances are then obtained by substi- 
tuting the corresponding values of the correlation 
coefficients Into Equation (15). 

The present study, however, uses Equations 
(16) and (18) during the entl'^e data reduction. 
The reason for this Is contained 1n the results of 
the sensitivity analysis presented in the next 
section. This analysis demonstrated that there Is 
not significant error magnification In the data 
reduction due to the correlation terms. 



4. Results of Hot-Wire Measurements 


Expressions for u'w' and v'w' can also be obtained 
in a similar manner and are given In Reference 22. 

1.2 Calculation of Covariances 


Dvorak and Syred^® used a DISA time correlator 
(55A06) to find the correlation coefficients 
between the velocity fluctuations In the three 
directions. King's approach Is to use the Infor- 
mation obtained by all six orientations and devise 
a mathematical procedure to calculate the covari- 
ances. 


Covariances are calculated using the relation- 
ship; 






(15) 


where y? 7 is the correlation coefficient between 
1 j 

the two cooling velocities Z^ and Zy By defini- 
tion, the absolute value of the correlation coef- 
ficient Y7 7 is always less than 1. 

King^® made certal' jssumptlons to calculate 
the covariances. However, he observed that at 
times the calculated value of the correlation 
coefficient ^s greater than one at which Instance 
he assigned previously fixed values to the corre- 
lation coefficients. He argued that If two wires 
are separated by an angle of 30 degrees, the fluc- 
jating signals from the wires at the two loca- 


The six-orientation hot-wire technique was 
employed to measure the turbulence quantities for 
swirling and non- swirling flow conditions In the 
confined jet facility described earlier. Also, an 
extensive sensitivity analysis of the data reduc- 
tion was conducted to assist the estimation of the 
uncertainties In the output quantities. 

4.1 Uncertainty Analysis 

The uncertainty analysis Includes a deter- 
mination of the sensitivity of the s1x-or1entat1on 
hot-wire data reduction to various Input param- 
eters which have major contributions In the 
response equations. Pitch and yaw factors (G and 
K) are used In the response equations described in 
Section 3 in order to account for the directional 
sensitivity of the single hot-wire probe. Figure 
6 shows the pitch and the yaw factors plotted 
against the hot-wire mean effective voltage. Both 
the pitch and yaw factors are functions of the 
hot-wire mean effective voltage, but the yaw 
Factor Is far more sensitive. A one percent 
Increase In the hot-wl.e voltage reduces the 
pitch factor by 1.3 percent and the yaw factor by 
56 percent. For the present study, the values of 
these factors are chosen at an average hot- 
wire voltage experienced In the flowfleld. 
TMs was a pprop rl ate since the output quantities 

(u, u' , u'v', etc) are only weakly dependent on 
rms 

the value of K. This can be seen In the data of 
Table 2 which summarizes a sensitivity analysis 
performed on the data reduction program at a 
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representative position In the flowfield. 

Table 2 demonstrates the percent change in the 
output quantities for a 1 percent change in most 
of the important input quantities. For the data 
presented in this table only quantities calculated 

from the probe orientation combination Z^, Zg and 
are used, for simplicity. In this swirling 
flow Zg was the minimum of the 6 mean effective 

cooling velocities. King^® has argued that the 
probe orientation combination approximately cen- 
tered around the minimum effective cooling velo- 
city produces more accurate estimates of calcu- 
lated turbulence quantities, than do the other 
orientation combinations. 


tions of data from adjacent wire orientations. 
One measure of the uncertainty in the output quan- 
tities can be obtained by examining the variance 
in these quantities calculated from the six dif- 
ferent position combinations. Table 3 shows these 
comparison data^for a representati ve position in 
the flowfield. For each of the output quanti- 
ties, an ensemble mean x is calculated together 

with art_ ensemble standard deviation o. The 

ratio a/x is a measure of the uncertainty in the 
output quantity. In this table, MR stands for 

'not resolved', a problem that occurs when the 

data reduction program attempts to take the square 
root of a negative quantity. In addition, quan- 
tities which are more than three standard devia- 
tions outside the mean are rejected as spurious 
calculations. 


Effect of Input Hr«mtt«r$ on Turbulpncp 
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\ CHWGE IN TIME-MEM AKO TUNMLEmCC QUMTITIES 

1 OAfllHFTCR 

: 

\ CHAMGi IM 
PWAf«TC8 

u 

V 

M 


D 


B 


i 

f 

j 1 


♦iS.JO 

*0.66 

*4.48 

•15.75 

• ?.06 

♦2.75 

♦6.0 

♦51.43 

♦11.94 

i - 

•1 


-Z.2\ 

♦11.44 

-6.50 

•2.4? 

*1Z.M 

♦4.0 

♦14.29 

♦ 7.46 


♦ 1 

-lO.SA 

-0.36 

•8.50 

-1.R8 

♦7.07 

•9.54 

•6.0 

•54.29 

-11.94 



• 0.27 

-0.06 

♦0.14 

♦1.6} 

*0.13 

♦0.39 

♦2.0 

♦2.66 

♦ 1.49 


•1 

*0.05 

0.0 

♦ 0.14 

0.0 

-0-13 

♦ 1.57 

0.0 

0.0 

♦1.49 

^6 


-0.1€ 

•0.18 

-0.14 

-0.6} 

•1.03 

•1.06 

-2,0 

-5.71 

g.O 



-1.02 

0.0 

•I.Ol 

•l.O 

0.0 

•C.98 

-2.0 

•2.S6 

-1.49 

< 

♦1 

*0.01 

-0.04 

•O.Ol 

♦0.01 

0.0 

• :.oi 

0.0 

0.0 

0.0 

'Vq 

M 

•0.05 

0.0 

•0.14 

-0.13 

-0.13 

-1.77 

0.0 

-2.86 

♦1.49 


♦ 1 

♦0.?I 

*0.01 

*0.05 

•1.63 

♦0.13 

-0,>9 

0.0 

•5.71 

♦ 1.49 


*1 

-0.16 

*0.1« 

-0.08 

*0.13 

0.0 

40.69 

•2.0 

♦2.86 

0.0 

.. .. 1 


It is not unusual in hot-wire anemometry to 
have the mean velocity components and turbulence 
quantities that are measured, be quite sensitive 
to changes in mean hot-wire voltage. For inter- 
precive purposes, the mean hot-wire voltage varia- 
tions can be thought of as being either errors in 
measuring the mean voltage, or shifts in the indi- 
vidual wire calibrations due to contamination or 
strain 'aging' of the wire. The data of Table 2 
demonstrate that the most serious Inaccuracies in 
the measurement and data reduction technique will 
be in the estimates of turbulent shear stresses, 

the most inaca.rate output term being u’w'. 

As already discussed in Section 3, an ad hoa 
assumption is made regarding the numerical values 
of the correlation coefficients used in the deduc- 
tion of time-mean and turbulence quantities. The 
results of the sensitivity analysis (Table 2) show 
the time-mean and turbulence quantities to be 
relatively insensitive to variations in the cor- 
relation coefficients. Therefore, the major ad 
';,v assumption made in the technique does not seem 
to have a great effect on the c/'r-t quantities 
compared to the effect of other input quantities. 

As mentioned earlier, turbulence quantities 
the output) can be calculated from six combina- 


T<8l.t I 


Scpttwr Awn) the Turbgtenc* Quantities When Solved by 
Si « different Cc^tnetions 


TUtAUUNCe 

TUR8ULENCC QJMriTf SOlVEO 87 SI< C0m81MAT10mS 

«AN 

S74NOMO 


QUMflTT 

1.2,3 

2.3.4 

1.4.5 

4.5.6 

5.6,1 

6.1.2 

■ 

OeVUTlOM 

a 

a/a 

"'“o 

0.21 

0.20 

0.21 

0.21 

0.L9 

0.16 

0.20 

0.01 

0.06 

v/u 

0 

0.10 

MB* 

0.11 

0.17 

0.17 

0.17 

0.14 

0.04 

0.26 

«/w 

0 

0.40 

0.39 

0.39 

0.3B 

0.37 

0.40 

0.39 

0.01 

0.03 


0.14 

C.I4 

0.14 

0.07 

0.08 

0.06 

O.U 

0.03 

0.31 


0.06 

0.11 

o.u 

0.08 

0.06 

0.09 

0.09 

0.02 

0.23 


0.13 

0.16 

0.10 

O.U 

0.10 

0,12 

0.12 

0.C2 

0.20 

77 t .> 

MR* 

MR' 

0.012 

MR* 

0.005 

0.004 

0.007 

0.004 

0.62 


0.002 

O.OlO 

0.002 

0.004 

0.004 

0.008 

O.OOS 

0.003 

0.72 


0.003 

MR* 

0.003 

0.003 

0.007 

O.OOl 

0,003 

0.002 

0.56 
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The data in Tables 2 and 3 can be used to pro- 
duce estimates in the uncertainties of the calcu- 
lated turbulence quantities. The data suggest 
that uncertainties on the order of 5 percent are 
to be expected in the mean velocity component 
estimates. Normal turbulent stress estimates 
(u'^i^g, etc.) have uncertainties on the order of 
20 to 30 percent and turbulent shear stress esti- 
mates are significantly hight , although most of 
this is a consequence of taking a product of terms 
such as u' and v' . 

These uncertainty estimates are considered to 
be somewhat conservative. More accurate estimates 
are quite difficult to obtain because, to our 
knowledge, similar measurements have not been per- 
formed with any other instrumentation system in 
this geometry flowfield. Also, comparisons of 
several representative points with independent 
measurements suggest that the ensemble averages 
estimates are typically in closer agreement than 
are selected sets of three orientations. There- 
fore, all turbulence estimates presented in this 
paper are calculated from ensemble averages of six 
groups of three adjacent wire orientations. Any 
data not resolved are not included in this avera- 
ging. This approach represents a departure from 
the technique developed by King*® who typically 
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x/D=1.5, r/0=0.25, ♦“SS deg (swirling flow). 




selected one group of three orientations from 
*«h1ch to calculate his turbulence estimates. 


4,? Results of Flowfleld Surveys 

Radial distributions of time-mean velocities, 
turbulent normal stresses and shear stresses are 
obtained for both nonswirling and swirling con- 
ditions, at various axial locations 1n the flow- 
field. 

Nonswirling Flow . In the confined jet, the 
experiments have Been conducted with expansion 
angles of 90 degrees (sudden expansion) and 45 
degrees (gradual expansion) and the results for 
both cases are presented In Reference 22 . In the 
Interest of brevity, only the data for a 90 degree 
expansion are presented here. 

Figure 7 shows the radial distribution of 
time-mean axial and radial velocity components at 
various axial locations. The axial velocity dis- 
tributions are compared with a similar study per- 
formed by Chaturvedl^ with a crossed hot-wire 
orobe. Because of the Inability of the slx-orlen- 
tation hot-wire technique to determine the sense 
of the flow direction In a nonswirling flow, the 
presence of the corner recirculation zone was 
observed by a sudden Increase in the axial velo- 
city closer to the wall . 

Figur, B shows the radial distribution of 
axial and radial components of the turbulence 
intensity at various axial locations In the con- 
fined jet flowfield. These turbulence intensity 
components are compared with Chaturvedi's measure- 
ments^ and reasonable agreement is found. In 
fact, the agreement In most cases is better than 
the uncertainty estimates derived from the data 
reduction sensitivity analysis. 

Included on Figure 8 are measured turbulent 
shear stress component (u'v'A)^) profiles for the 

nonswirling confined jet. For the most part, 
these measurements are In reasonable agreement 
with those made by Chaturvedl^ with a crossed wire 
probe. The two significant exceptions to the good 
agreement occur at the furthest upstream and fur- 
thest downstream locations. Upstream, at x/D = 
0.5 the shear layer Is very thin and, therefore, 
matching data from several wire orientations 
obtained at somewhat different times may be prac- 
tically difficult. We believe the overly large 
measured turbulent shear stress on the centerline 
at the furthest downstream station (x/D = 3.0) to 
be a consequence of the transient nature of the 
flow. The recirculation regions In the confined 
jet oscillate somewhat at a low frequency, likely 
characteristic of the main acoustic modes In the 
tube. These large scale oscillations can have 
significant correlated velocity fluctuations (such 

as 7V^). 

Swirling Flow The measurements performed In 
the swi r1 1 ng fl ow are with a = 90®, 4i = 38®, and 
x/D = 0.5, 1.0, and 1.5. The object of these 
limited number of experiments was to evaluate the 
reliability and accuracy of the six-orientation 
hot-wire technique before making extensive use of 
the technique. 
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The hot-wire results In the case of time-mean 
axial and azimuthal (swirl) velocities shown In 
Figures 9 and 10, are compared with five-hole 
pitot probe measurements performed by Rhode^ In 
the same experimental facility. Agreement among 
the two studies Is fairly good. King*^® suggested 
a method to determine the sense of the axial velo- 
city. He advised comparing the magnitudes of Z3 
and Zg given by Equations 4. In the present flow- 
field, the swirl velocity Is always positive and 
the two equations giving Z3 and Zg differ only In 
the sign of axial velocity. Therefore, when Zg 
Is greater than Z3, the axial velocity 1s nega- 
tive, otherwise It Is positive. With the proper 
sense being assigned to the x velocity mean com- 
ponent, the presence of central torlodal recir- 
culation zone Is evident in the results of both 
measurement techniques. 

Figure 11 shows the radial distribution of the 
time-mean radial velocity at various axial loca- 
tions for a swirl vane angle of 38 degrees and 
wall expansion angle of 90 degrees. Flow visuali- 
zation and five-hole pitot probe measurements per- 
formed in Rhode's study® show the time-mean radial 
velocity to be negative at axial locations greater 
than x/0 = 0.5. In spite of the Inability of the 
six-orientation hot-wWe technique to determine 
the sense of the radial velocity, the data are 
presented with the appropriate sign change. There 

Is a reasonable agreement among the two studies In 
measurements of time mean radial velocities except 
at the initial measurement station. 


Figure 12 shows the radial distribution of 
axial, radial and azimuthal turbulent Intensities 
at three axial locations presented. At axial 
locations closer to the Inlet of the confined jet, 
the axial turbulence Intensity Is fairly high, up 
to 32 percent for x/D = 0.5 which Is due to the 
large axial velocity gradients closer to the wall. 
However, in the case of radial turbulence inten- 
sity, the profiles are rather flat. The mean azi- 
muthal velocity also experiences sudden changes in 
gradients and, hence, the outcome Is a large azi- 
muthal turbulence Intensity closer to the wall at 
x/D = 0.5. 


Figure 13 shows the shear stresses u v /U„, 

2 2 0 

u'w’/Uq, andv'w'/U^as a function of radial and 

axial distance. The sensitivity analysis showed 
that we should expect large uncertainties asso- 
ciated with evaluation of turbulent shear stresses 
using the six-orientation technique. Therefore, 
the reliability of the profiles of these shear 
stresses shown in Figure 13 Is uncertain at 


this time. Nevertheless, stresses u v /U^ and 

— 2 0 

u'w'/U^ are found to have large values closer to 

the wall that one would expect due to steep axial 
and azimuthal velocity gradients. The fact that 
we have found no other measurements of this type 
In a swirling, recirculating flow attests to the 
fact that accurate measurements In such a flow are 
quite difficult. 


Cl osure 


The six-orientation hot-wire technique Is a 
relatively new method to measure time-mean velo- 
city components and turbulence quantities In com- 
plex three-dimensional flowflelds. Applied In 
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Figure 7. Radial Distributions of Time-Mean Axial 
and Radial Velocity Components In the 
Non-Swirling Confined Jet; (Note the 
Difference In Velocity Scales) 
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Figure 8. Radial Distributions of 

Turbulent Intensities u' /u^, v' /u„ 
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Figure 9. Radial Distribution of Time-Mean ^lal 
Velocity In Swirling Confined Jet for 
a ■ 90* and 9 ■ 38’ 
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Figure 10. Radial Distribution of Time-Mean 
Azimuthal Velocity In Swirling 
Confined Jet for a • 90* and 9 - 38* 


this study to nonreacting a x1 symmetric flowflelds, 
measurements of time-mean and root-mean-square 
voltages at six different orientations contain 
enough information to obtain the time-mean velo- 
cities, turbulence intensities and shear stresses. 
At each location In the flow, there are six dif- 
ferent values of each of the above quantities that 
can be obtained using six sets of measurements of 
three adjacent orientations. Ensemble averages of 
the output quantities from the six combinations of 
data appear to produce estimates with the best 
agreement with independent measurements. 


Flowfleld surveys of both swirling and non- 
swirling confined Jets have been made with the 
si x-or1entatlon single hot-wire technique. These 
measurements have been used to calculate estimates 
of the mean velocity components and the normal and 
shear turbulent stresses. Where Independent data 
exist, comparisons have been made which demon- 
strate the reliability of the technique. 

In addition, a sensitivity analysis of the 
data reduction technique has been conducted wMch 
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forms the major Ingredient In the uncertainty 
analysis. It Is demonstrated that the largest 
uncertainties are to be expected in the turbulent 
shear stress estimates. Nevertheless, in non- 
xwirlinq flows the measured shear stresses are in 
reasonable agreement with previous measurement! 
made with a crossed-wire probe. In swirling flow, 
previous similar measurements have not been found. 
Consequently, the universal accuracy of the mea- 


sured technique applied to swirl flows is still an 
open question. 
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Figure 11. Radial Distribution of Time-Mean 


Radial Velocity in Swirling Confined 
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Fi^re 12. Radial Distribution of Turbulent Shear Stresses u'v'/u*. u'w' /u*, 
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